Photonic realization of the quantum Rabi model 
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We realize a photonic analog simulator of the quantum Rabi model, based on light transport 
in femtosecond-laser- written waveguide superlattices, which provides an experimentally accessible 
testbed to explore the physics of light-matter interaction in the deep strong coupling regime. Our 
optical setting enables to visualize dynamical regimes not yet accessible in cavity or circuit quantum 
electrodynamics, such as bouncing of photon number wave packets in parity chains of Hilbert space. 
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The quantum Rabi model [TJ [2] is a milestone in the 
history of quantum physics. It describes the simplest 
interaction between quantum light and matter, i.e. a 
two-level atom coupled to a single quantized electromag- 
netic mode, and applies to a great variety of physical 
systems, including cavity and circuit quantum electro- 
dynamics (QED), quantum dots, polaronic physics and 
trapped ions. In spite of its simplicity, the integrabil- 
ity of this model has been proven only recently [3 , and 
the experimental exploration of some extreme dynami- 
cal regimes of the Rabi model is still missing. In most 
experimental conditions, such as in cavity QED, simul- 
taneous creation or annihilation of an excitation in both 
atom and cavity mode is extremely unlikely, and the co- 
herent atom-field dynamics is described by the Jaynes- 
Cummings (JC) model [4], which is obtained from the 
quantum Rabi model via the rotating-wave approxima- 
tion (RWA). The JC model has been very successful to 
accurately predict several experimental phenomena, such 
as the observation of the vacuum Rabi mode splitting 
in Alkali atoms [5] and vacuum Rabi oscillations in Ry- 
dberg atoms [6 . Recent achievements in circuit QED 
enabled to explore the ultrastrong coupling regime of 
light-matter interaction [7H9] , in which effects of counter- 
rotating terms can not be neglected, but can still be 
taken into account by means of a perturbative approach. 
This has led to the observation of the Bloch-Siegert shift, 
i.e. a noticeable shift in the resonance frequency of the 
driven atom [9]. Recently, an even stronger coupling in 
the field-atom interaction, the so-called deep strong cou- 
pling (DSC) regime, has been theoretically investigated 
[TO] . This regime includes the highly counterintuitive si- 
multaneous excitation or de-excitation of both the atom 
and the field. Concepts like Rabi oscillations should be 
abandoned, and new extreme physical phenomena, such 
as bouncing of photon number wave packets along par- 
ity chains in Hilbert space, come into play [TO]. Although 
cutting-edge experiments that can access this novel light- 
matter coupling regime are desirable, the current status 
of cavity or circuit QED does not allow to reach the DSC 
regime yet. Creating tabletop analog systems to simulate 



in the lab extreme dynamical regimes or phenomena in 
the matter, which are not accessible in an experiment, 
has attracted an increasing interest in recent years (see, 
for instance, [TTUT3] and references therein). The possi- 
bility to simulate the DSC regime of light matter interac- 
tion, using either classical or quantum systems, has been 
recently proposed in Refs.p~4l [15]. In particular, in [14] 
it was shown that propagation of classical light beams in 
curved photonic superlattices provides a classical realiza- 
tion in Hilbert space of the quantum Rabi model which 
can access all kinds of light-matter coupling regimes. 
In this Letter we realize a classical simulator of 
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FIG. 1. (Color online) (a) Schematic of a qubit interact- 
ing with a single bosonic mode, (b) Photonic realization of 
the quantum Rabi model in Hilbert space based on a semi- 
infinite curved binary waveguide array, (c) Variation A/3 in 
the propagation constants f3 of the waveguides, induced by 
the axis curvature: actual trend without compensation due 
to a non-uniform waveguide spacing (diamonds) and target 
linear trend achieved after writing parameters compensation 
(solid line), (d) Top-view microscope image of the curved 
waveguide array. 
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the quantum Rabi model, based on light transport in 
femtosecond-laser-written waveguide superlattices, and 
provide experimental evidence of the physics of light- 
matter interaction in the DSC regime. 

The Hamiltonian of the Rabi model, describing the 
interaction of a two-level atom (or a qubit) with a quan- 
tized field (or oscillator) (FigjT^i), reads [2j|3] 

H = (hjjj$/2)cF z + houa^a + Hg(&+ + a_)(a + a^), (1) 

where a and a) are the annihilation and creation opera- 
tors of the quantized oscillator with frequency lj, a z = 
\e)(e\-\g)(g\, cr + = \e)(g\ and <r_ = \g)(e\ are Pauli oper- 
ators associated to the two-level atom with ground state 
excited state |e), and transition frequency ujo, and g 
is the coupling strength. For a weak coupling g <C uJo and 
near-resonant excitation uj ~ u;o, the anti-resonant terms 
a + a) and <r_a entering in H can be neglected (RWA), 
and the Rabi model reduces to the JC model of quan- 
tum optics [2|. From a physical viewpoint, neglecting the 
anti-resonant terms means that processes corresponding 
to simultaneous excitation or de-excitation of the atom 
and of the quantized field are unlikely. In this regime, the 
dynamics in Hilbert space, spanned by the states \g)\n) 
and \e)\n) describing n quanta of the field with the atom 
in the ground or in the excited state, splits into an infinite 
sequence of state doublets, and phenomena like vacuum 
Rabi oscillations or collapse and revival of populations 
can be observed. Without the RWA, the dynamics in the 
Hilbert space is more complex and splits into two parity 
chains [10]. In fact, let us expand the state vector \i/)(t)) 
of the atom-field system on the {\g)\n}, \e)\n)} basis as 
I^CO) = E^LoKWI e )l n ) + & nWb)|^)], and let us intro- 
duce the amplitudes c n (t) and f n (i), defined by c n (t) = 
Q>n •> fn(t) = b n for n even, and c n (t) = b n , f n (t) = a n for 
n odd. Then, the temporal evolution of the amplitudes 
c n and f n decouples into two parity chains as follows [TO] 

dc n , . ujq , . 

i—jjT = ^nC n +l + Kn-lCn-1 + (-1) —C n +nUJC n (2) 
= ^n/n+l + Kn-lfn-1 (-1)" y/n + nuf n (3) 

(n = 0, 1, 2, 3, ...), where we have set n n = g\fn + 1. 

As noticed in [14], Eqs.(2) and (3) are analogous to 
the coupled-mode equations describing light transport 
in two decoupled semi-infinite binary photonic lattices 
in the tight-binding approximation with a superimposed 
transverse index gradient and with non-uniform coupling 
constant n n between adjacent waveguides n and (n + 1). 
In the optical analogue, the temporal evolution of the 
amplitude probabilities c n and f n is mapped into spa- 
tial propagation of the light amplitudes in the various 
waveguides (see Fig. lb). The transition frequency uo of 
the two-level atom corresponds to the propagation con- 
stant mismatch of alternating waveguides of the super- 
lattice, whereas the frequency uj of the quantized mode 
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FIG. 2. (Color online) (a) Measured and (b) predicted light 
intensity distributions in the waveguide array for uo — 0, 
showing the periodic bouncing of photon number statistics 
P(n,t) in Hilbert space, (c) Evolution of population P e of 
level \e) (upper panel) and of the revival probability P r (lower 
panel) (crosses: experimental points; solid line: theoretical 
model). 



is defined by the waveguide propagation constant gra- 
dient superimposed to the superlattice. Note that the 
two parity chains of Eqs.(2) and (3) are simply obtained 
from each other by changing the sign of ujq. The struc- 
ture of the optical superlattice that realizes one of the 
two parity chains is shown in Fig{T]3. It comprises a 
semi-infinite planar array of curved waveguides with en- 
gineered distances d n in order to implement the desired 
coupling coefficients n n . The circular curvature of the 
waveguides introduces a transverse refractive index gra- 
dient leading to uj = 2i\n^ yU d n / \RX) , where n e ^^ n is 
the effective refractive index of the n-th waveguide, A is 
the wavelength of light and R the radius of curvature. 
Since the distances d n between different waveguides are 
not uniform in the array (Figjlj}), the value of uj also 
varies with the waveguide number n (see diamonds in 
FigfTJ^). To achieve a uniform uj in the array a com- 
pensation is required in the effective index n e ^ y7l of the 
waveguides, yielding a uniform n e ^^ n d n product. In ad- 
dition to this uniform transverse index gradient, the ujq 
detuning is implemented by a further modulation of the 
effective refractive index in alternating waveguides. The 
present design of the photonic device differs from that 
proposed in [14], where a three-dimensional layout was 
considered. In an actual implementation that geometry 
suffered from second-neighbours interactions and difficul- 
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FIG. 3. (Color online) (a) Same as Fig. 2c, but for a waveguide array with cjo= 0.04 mm -1 , (b) Early evolution of revival 
probability P r and populations P e (for cjo > 0) and P g — 1 — P e (for cjo < 0) as cjo is tuned from negative (initial state 
|^(0)) = |#)|0)) to positive (initial state |^(0)) = |e)|0)) values. Solid curves: theory; crosses: experimental points. The dashed 
horizontal lines in the plots show, as reference levels, the minimum values of revival probability and populations attained in 
the first bounce interval for the ujq — case. 



ties in simultaneous imaging of waveguides at different 
depths. 

The photonic structure shown in FigJT}} is manufac- 
tured by femtosecond laser waveguide writing on a fused 
silica substrate. The second harmonic of an Yb-based 
femtosecond laser (FemtoREGEN, HighQLaser GmbH), 
delivering 400 fs pulses, is used for the writing process. 
An optimal processing window was found at 20 kHz repe- 
tition rate, 300 nJ pulse energy and 10 mm/s translation 
speed. The laser beam is focused 170 jam below the glass 
surface by a 0.45 NA, 20 x objective. Figure [T]i shows 
a top-view microscope image of a typical manufactured 
array of 15 curved waveguides (R = 650 mm). Note 
that the separation between the waveguides is not uni- 
form, ranging from 6.6 /im to 14 /am. The modulation of 
the effective refractive index n e ff jn , requested for achiev- 
ing both a uniform index gradient uj and the alternating 
mismatch co>o, is obtained by exploiting the refractive in- 
dex change dependence of the guiding core on the writing 
speed. Calibration of this dependence provided a rate of 
dn e ff^ n /dv = 1.5 xlO -5 s mm -1 , thus requiring a varia- 
tion of the writing speed in the range 10-14 mm s _1 to 
achieve a uniform uj and a further modulation up to ±0.3 



mm s 1 to implement the detuning ujo. 

Femtosecond laser writing in fused silica creates color 
centers that provide fluorescent emission at about 650 nm 
when light at 633 nm is propagated in the waveguide [16] . 
Top-view imaging of the fluorescence signal is a conve- 
nient way to visualize and quantitatively estimate the 
light distribution along the waveguide array [17], reject- 
ing the background light by a notch filter at 633 nm. 

In all the experiments presented below an array of 15 
curved waveguides was manufactured with parameters 
uj = 0.23 mm -1 and g = 0.15 mm -1 . The condition 
g — uj breaks the RWA and allows the investigation 
of the DSC regime. A first experiment addressed the 
limiting case ujo = 0, which corresponds to degenerate 
qubit levels. In this case, H can be diagonalized by a 
Lang-Firsov transformation [18], and the energies turn 
out to be equally spaced by cj, yielding a strictly peri- 
odic temporal dynamics of photon number distribution 
P(n,i) and populations with period T = 2tt/uj. Figure 
[2^i shows the top- view image of light distribution in the 
array, which provides the evolution of P(n, t). In the 
experiment, the n = waveguide is excited at the in- 
put plane (left side in Fig. 2a), corresponding to the 
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system in the initial state |^(0)) = |e)|0). A periodic 
bouncing of the photon distribution can be clearly ob- 
served, in very good agreement with the expected one 
(Fig {2^). The qubit population of the |e) level is given 
by P e (t) = S^Lo l c 2n| 2 , while the revival probability 
is defined as P r (t) = |(^(t)|^(0))| 2 = |c | 2 . The evo- 
lution of these two quantities can be readily retrieved 
from the light intensity map of Fig(2^i, and is depicted 
in Figj2]3. Note that the measured behaviors of popula- 
tion and revival probability closely follow the theoretical 
curves, with a period T ~ 27.3 mm. 
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FIG. 4. (Color online) (a) Measured light intensity distribu- 
tion in a waveguide array for cjo = -0.08 mm -1 (left panel) 
and cjo = +0.08 mm -1 (right panel), (b) Evolution of the 
mean photon number (n) for cjo = 0, ±0.08 mm -1 . 

When ujq 7^ 0, the energies of H are not equally spaced, 
and the bouncing dynamics ceases to be periodic [TO] . 
However, for uo smaller than cj, an imperfect bouncing 
dynamics can be still observed in the early stage of the 
evolution. As an example, Fig(3^i shows the evolution of 
P e and P r observed in an array with ujo = 0.04 mm -1 , 
together with the corresponding theoretical curves, for 
the system in the initial state |^(0)) = |e)|0). Smearing 
out of the periodic dynamics at successive bounces can 
be appreciated in this case. For ujq ^ 0, the evolution 
of both revival probability and populations, besides los- 
ing the exact periodic behavior, is sensitive to the sign of 
co>o- In fact, according to Eqs.(2) and (3) the sign of Uq 
defines two distinct initial states of the qubit: the qubit 
is in the excited state for c^o > 0, whereas it is in the 
ground state for ujq < 0. Such different initial conditions 
excite two different parity chains in Hilbert space, leading 
to different field-atom dynamics for the same initial pho- 
ton number distribution. This is clearly shown in Figj3fD, 
which depicts the early evolution of both revival proba- 
bility P r and populations P e (for ujq > 0) and P g = 1 — P e 
(for ujq < 0) as ujo is tuned from negative to positive val- 
ues, with no field quanta in the initial state. Note that, 



as cuo is increased from negative to positive values, the 
minimum value of both the population and revival prob- 
ability, reached within the first bounce, decreases. Such a 
behavior can be physically understood by observing that, 
when the qubit is initially in the \g) state (i.e. for uoq < 0) 
and because of the absence of field quanta in the initial 
state, the decoupled atom-field system is in its lowest en- 
ergy state and excitation of field quanta and atom is thus 
less probable for an assigned coupling g with respect to 
co>o > 0- The possibility to create excitation in the atom 
and field starting from the initial state \ip(0) = |#)|0), as 
demonstrated in the left panels of Fig. 3b corresponding 
to uo < 0, is a rather amazing prediction of the DSC 
regime of the Rabi model. Such a situation corresponds 
to creation of photons from an atom initially in its ground 
state which is strongly interacting with the vacuum state 
of the field. In our ordinary wisdom of atom-field inter- 
action, such a possibility would be obviously prevented 
because of energy conservation. In the DSC regime we 
actually do not violate energy conservation, because the 
energies of excitations (Hljo and two) are of the same order 
of magnitude or smaller than the "coupling" energy (%). 
To better highlight the different dynamics of the cou- 
pled field-atom system corresponding to the initial states 
|^(0) = \g)\0) and |^(0) = |e)|0), in Figgwe show the 
detailed evolution of photon number distribution P(n, t) 
for ujo — -0.08 mm -1 and ujo = 0.08 mm -1 , together with 
the average number of created photons (n). Note that in 
both cases (n) follows the bouncing dynamics of P(n,i). 
The maximum value of (n) within the bouncing cycle is 
noticeable larger in the ujq > case. This behavior is 
basically due to the fact that in this case the atom is 
initially in the excited state and such an extra energy is 
available to create photons. 

In conclusion, we have experimentally investigated the 
physics of light-matter interaction in the deep strong cou- 
pling regime using a photonic analog simulator of the 
quantum Rabi model based on light transport in engi- 
neered waveguide superlattices realized by femtosecond 
laser writing. Our optical set up provides a simple tool 
to visualize in the lab extreme dynamical regimes of light- 
matter interaction not yet accessible in cavity or circuit 
QED systems, and could be further extended to simulate 
non-Hermitian extensions of the quantum Rabi model or 
related models [19] by introduction of optical gain and 
absorption in the guiding structure. 
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